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^ ■ Two examples of the use of differential geometry in plasma physics are given: The first is 

Q , the computation and solution of the constraint equations obtained from the Riemann met- 

^ ric isometry of the twisted flux tube. In this case a constraint between the Frenet torsion 

Q ■ and curvature is obtained for inhomogeneous helical magnetic flux tube axis. In the second 

c/3 . one, geometrical and topological constraints on the current-carrying solar loops are obtained 
by assuming that the plasma filament is stable. This is analogous to early computations 

^ ' by Liley [(Plasma Physics (1964)] in the case of hydromagnetic equilibria of magnetic sur- 

H : 

. . . faces. It is shown that exists a relationship between the ratio of the current components 

along and cross the plasma filament and the Frenet torsion and curvature. The computations 
are performed for the helical plasma filaments where torsion and curvature are proportional. 
The constraints imposed on the electric currents by the energy stability condition are used 
to solve the remaining magnetohydrodynamical (MHD) equations which in turn allows us to 
compute magnetic helicity and from them the twist and writhe topological numbers. Mag- 
netic energy is also computed from the solutions of MHD equations. PACS numbers: 

02.40.Hw-Riemannian geometries 
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I Introduction 



The tools of differential geometry and ,in particular, Riemannian geometry [1] have been 
proved very useful in handling problems in areas of physics ranging from plasma astrophysics 
[2] to Einsteins general relativity [3]. This common feature can be better appreciated when 
one notices that in one of the problems we shall address in this paper, the isometry which is 
the vanishing of the variation of the metric Qik , or 6gik = = 1,2,3) either in gravita- 
tional spacetime, or as we shall see here, in the magnetic flux tube axis geometry. Another 
application of differential geometry and topology to plasma physics is the investigation of 
topology and geometrical constraints in fast kinematic dynamos [4, 5] where writhe number 
and twist were computed from magnetic helicity, or as shall address in this paper, to the inves- 
tigation of the stability of solar loops in plasma astrophysics. Riemannian geometry has also 
been used with great success to investigate the tokamaks and stellarators metric oscillations 
and the instabilities of confined plasmas [6]. In this paper we show that when one imposes 
the constraint of twisted magnetic flux tube [7] metric isometry a relation between Frenet 
curvature [8] and inhomogeneous helical solar loops are obtained. The helicoidal form of the 
device implies that Frenet torsion of magnetic field lines is compulsory. Earlier Lileys [9] has 
computed the current conditions for static hydromagnetic equilibria. In plasmas the perturba- 
tion of energy 5W of magnetic surfaces is computed, and one is able to say if the magnetic 
surface is stable or unstable according the 6W > or 5W < respectively. The extreme 
condition 5W = is considered to be the case of stable equilibria in the case understand 
the mechanism which allow the highly twisted coronal magnetic flux tube to emerge from the 
solar surface and produce those beautiful solar flares and loops. Recently Berger and Prior [10] 
have presented a detailed discussion of twist and writhe of open and closed curves presenting 
a formula for the twist of magnetic filaments in terms of parallel electric currents, s paper we 
also compute the writhe and twist numbers investigated previously by Moffat and Ricca [11] 
and Berger and Field [12] to obtain an expression for the writhe number for the magnetically 
perturbed vortex filamentary twisted structure. To obtain the writhe number we make use of 
the helicity local expression for the magnetic field proportional to the magnetic field itself and 
decompose the magnetic vector field along the Serret-Frenet frame in 3D dimensions. This 
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allows us to solve scalar MHD equations to obtain constraints on the global helicity expression 
which will allow us to connpute the writhe number. The knowledge of tilt , twist and writhe of 
solar filannents for exannple has recently helped solar physicists [13] to work out data obtained 
from the vector magnetograms placed in solar satellites. This is already an strong motivation 
to go on investigating topological properties of these filamentary twisted magnetic structures. 
Earlier the Yokkoh solar mission has shown that the sigmoids which are nonplanar solar fila- 
ments are obtained due to the action of electric currents along these filaments which was used 
recently [14] as motivation to investigate current-carrying torsioned twisted magnetic curves. 
The paper is organized as follows: In section 2 we compute the conditions imposed on the 
current-carrying solar loops by their stability. Section 3 addresses the computation of writhe 
and twist numbers from the magnetic helicity while in section 4 we investigate the Riemann 
isometry constraint on the twisted thin filament and solve the scalar equations obtained for 
the curvature and torsion of the magnetic axis of the flux tube. In section 5 we present the 
conclusions. 



II Stability of current-carrying solar loops 

Let us now start by considering the MHD field equations 

V.B = (II.l) 

V.J=0 (II.2) 

where a is the magnetic twist and the magnetic field B along the filament is defined by the 

expression B = Bgt and j = jst + inn + jj) is the electric current density where js and are 
respectively the components along and cross to the solar filament all along its extension, and 
Bg is the component along the arc length s of the filament. The vectors t and n along with 

— * 

binormal vector h together form the Frenet frame which obeys the Frenet-Serret equations 



F = Kfi (II.3) 

n' = -Kt + rb (II.4) 

b' = -rn (II.5) 
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the dash represents the ordinary derivation with respect to coordinate s, and K{s,t) is the 
curvature of the curve where k — R~^. Here r represents the Frenet torsion. We follow the 
assumption that the Frenet frame may depend on other degrees of freedom such as that the 
gradient operator becomes 

d d 

The other equations for the other legs of the Frenet frame are 

d ^ 

—t^0nsn+[^b + r]b (II.7) 
d ^ ^ ^ 

—n = -Onst - (divb)b (II.8) 
on 

= + T]t- (divb)n (II.9) 

on 

^t^0j-[Qn + r]n (11.10) 

d 

fi^\fl^ + r]t- K + (divn)b (11.11) 

ob 

—b = -Obst -[k+ (divn)]n (11.12) 
do 

Now that we have this mathematical machinery at our disposal, let us consider the Bernstein 
et al [15] plasma energy stability relation 

5W=^J ds[{Q + n.tjxnf + -fp{dive)^ - 2{n.e)'^jxn.BVn)] (11.13) 

where 

Q^curl{exB) (11.14) 
and 7 is the ratio of the specific heats. According to Lileys [14] if one puts 

S ^jxn.{BVn) (11.15) 

and if there if there is no discontinuities in the equilibrium values, then a sufficient condition 
for stability , i.e, 5W>0, is 

S<0 (11.16) 
substitution of the expressions for j and for B 

jxn^jbt-j,b (11.17) 



into this expression yields the following geometrical constraint in the solar loop 

[Mti{s)+T{s)js]B>Q (11.18) 

By assunning that B > one obtains in the stable equilibrium case SW = or 5" = which 
reads 

^ = -- (11.19) 

js K 

Note that in planar (torsionless) solar loops the normal current jb vanishes. In helical solar 

loops where the ratio between torsion and curvature is constant or (3 = — ^ one has the 
following relation 

^ = -Z?!^ (11.20) 
as as 

Physically this means that as observe the solar loop along its extension the toroidal current 
increases while the normal current decreases as going to a focal point. A long and straight- 
forward computation allows us to write down the following expression 

V.J=0 (11.21) 

implies 

dsjs + [Obs + divb]js = -jb{divb + 9bs) (11.22) 
substitution of expression (11.19) into (11.22) allows us to obtain a differential equation for the 



current jg as 



dsjs + ebs[l--]js = (11.23) 



which yields the solution 

Js = Js'expl- f - -)ds] (11.24) 

Note that the current along the solar loop decays as the coordinate s— >oo. From this last 
expression one notes that 

js _ exp[- J Obsds] 

where the RHS inequality is valid since the initial current is greater than the resultant current. 
This inequality is equivalent to 

exp[— j Obsds] < exp[— j —ds\ (11.26) 



which yields the following constraint between the ratio between torsion and curvature and the 
factor 9bs as 

ebs > - (11.27) 

in all points of the solar loops which are possibly broken in the inflexionary points of the loops 
[7]. our computations we neglect divb and also consider that 9f,s varies very slowly along the 
solar loop. Now let us compute the magnetic field from the equation V.B — which becomes 

dsBs + [9bs + ens]Bs = (11.28) 
A simple solution of this equation can be obtained as 

B, = Boexp[ J {Obs + 9ns)ds\ (11.29) 

where Bq is an integration constant. Substitution into equation (11.17) assuming that the flow 
is geodesic along the filament and abnormality relation = along with 9i,s — 9ns — Done 
obtains 

S„ = BQ[{dwh)-^ - j Kds] (11.30) 
In the next section we shall compute the magnetic energy in terms of the magnetic vector 

— * 

potential A, which allows us to compute the helicity integral and the writhe number and twist 
of the solar loops. 



Ill Magnetic energy, Twist and Writhe of Solar Loops 

In this section we compute the magnetic energy , twist and writhe of the solar filament based 
on the vector potential, assuming that it obeys the Coulomb gauge V .A = and the definition 
B — VxA. The Coulomb gauge becomes 

dsAs + [Obs + 9ns]As = (III.31) 

assuming that A — A{s,n) together with the equations for the definition of B, taking A — 

As^nt y'\e\ds 

dnA + kA^O (III.32) 
[2r + Q„ + QsjA = (III.33) 
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which upon separation of variables A = $(s)6(n) one is able to find 

A{s, n) = Aoexp[-Kon + j {9bs + Ons)ds\ (III.34) 
Let us now compute the magnetic energy 

Eb^ — [ B^dV (III.35) 

which yields 

Eb = ^^^exp[-2{ebs + ens)s] (III.36) 

o 

Note that from above constraints — ^ and allowing 9ns to vanish one obtains the energy 

Eb = ^exp[-2(^).] (111.37) 

which is the energy of solar loops based on the assumption of its stability. Now let us compute 
the global helicity 

J A.BdV (III.38) 

which is equivalent to 

H ^ J [A(s,n).Bs]dV (III.39) 
substitution of the relations above yields 

H = -Tra^ J [2t + + QjA^rfs (III.40) 

Besides if one considers that the plasma flow in solar loops is geodesic the abnormality factor 
Qg vanishes and the magnetic helicity reduces to 

H = -na^ J [2t + (III.41) 

If one assumes that the magnetic vector potential varies also quite slow along the solar filament 
one has 

H = -na^A^ J [2t + Q„](is (III.42) 

The first integral on the RHS in (111.42) is the total torsion and considering that the helicity 
is proportional to the sum of twist and writhe and yet that the twist is proportional to the 
total torsion , we conclude that the writhe number Wr is proportional to the integral of the 
abnormality Q^. 



IV Riemann isometry of magnetic flux tube metric 

Since in the limit of thin magnetic flux tube one could represent it as a mathematical model 
for solar loop, it is interesting to argue if there is a similar way to constraint the magnetic flux 
tube as was done for the current-carrying solar loops in previous sections. A simple way to 
address and respond to this question is to resource to a simple tool of Riemannian geometry 
[7, 16] used in general relativity, which is called metric isometry and is represented by the 
expression 

8gik = (IV.43) 
If one applies the isometry hypothesis to the flux tube metric 

^ ^ ^2^^^ ^ K'^{s)ds'^ (IV.44) 

one obtains that 5grr — and that 

Sgee = [^9rr]5r = 2r5r = (IV.45) 
which yields r = and means we are close to the magnetic axis of the tube, while 

^9ss = [-Q:^9ss]Sr + [-Q^9ss]Ss (IV.46) 

which yields 

d d 

Sgss = [r{K{s)sine—e - —K{cose))]5s (IV.47) 

where we assume that 5r — Q or that we are over the magnetic axis of the flux tube. Here 
we also consider as in Riccas paper [8] that 9 — 9r — J rds. Applying the isometry condition 
above for this component one finally obtains the constraint equation for the Frenet curvature 
as 

= -To's + Co (IV.48) 

which after a simple integration yields 

d 1 

—K = --To's' + Cos + Ci (IV.49) 
OS 2 

This equation was obtained under the assumption that the magnetic axis of the twisted mag- 
netic flux tube torsion is constant t — To and that Co and Ci are integration constants. 



V Conclusions 



In conclusion, topological and geometrical constraints on solar loops were obtained through 
the use of the tools of differential geonnetry. Other kinds of isonnetries such as the Riennann 
curvature tensor Rij^i isometry SRijki — could be calculated for tokannaks and stellarators 
as well as for solar loops. This may appear elsewhere. We also hope that the analytical models 
discussed here may be useful to solar physicists in the checking their observational results in 
magnetic twist and writhe in solar activity regions. 
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